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Abstract. In this paper we prove a global existence theorem, in the 
direction of cosmological expansion, for sufficiently small perturbations 
of a family of spatially compact variants of the k = — 1 Friedmann- 
Robertson-Walker vacuum spacetime. We use a special gauge defined 
by constant mean curvature slicing and a spatial harmonic coordinate 
condition, and develop energy estimates through the use of the Bel- 
Robinson energy and its higher order generalizations. In addition to the 
smallness condition on the data, we need a topological constraint on the 
spatial manifold to exclude the possibility of a non-trivial moduli space 
of flat spacetime perturbations, since the latter could not be controlled 
by curvature-based energies such as those of Bel-Robinson type. Our 
results also demonstrate causal geodesic completeness of the perturbed 
spacetimes (in the expanding direction) and establish precise rates of 
decay towards the background solution which serves as an attractor 
asymptotically. 



1. Introduction 

In this paper we establish global existence and asymptotic behavior, in 
the cosmologically expanding direction, for a family of spatially compact, 
vacuum solutions to the 3+1 dimensional Einstein equations for sufficiently 
small perturbations of certain known "background" solutions. The back- 
grounds we consider are the spatially compactified variants of the familiar 
vacuum k = — 1 Friedmann-Robertson- Walker (FRW) solution, which exist 
on any 4-manifold M of the form (0, oo) x M, where M is a compact hy- 
perbolic 3-manifold (i.e., a manifold admitting a Riemannian metric with 
constant negative sectional curvature). 

Let 7 be the standard hyperbolic metric with sectional curvature —1 on 
M. Then (M, 7) given by 

M = (0,oo)xM 
7 = —dp <g> dp + /? 2 7 
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is a flat spacetime, locally isometric to the k = — 1 vacuum FRW model, 
which we shall call a hyperbolic cone spacetime. It has a big bang singular- 
ity as p \ but expands to infinite volume as p /* oo. The vector field pj^ 
is a timelike homothetic Killing field on (M, 7) so that these backgrounds 
are continuously self-similar. We shall be considering sufficiently small per- 
turbations of such hyperbolic cone spacetimes to the future of an arbitrary 
p =constant Cauchy surface under the additional topological restriction that 
(M, 7) be "rigid" in a sense that we shall define more fully below. The rigid- 
ity assumption serves to eliminate the possibility of making non-trivial but 
still flat perturbations of the chosen backgrounds. 

Our main result treats the vacuum Einstein equations on M and proves 
global existence in the expanding direction for initial data sufficiently close 
to data for (M, 7). More precisely, we show that the maximal globally hyper- 
bolic future vacuum development (M, g) of such data is causally geodesically 
complete and globally foliated by constant mean curvature (CMC) hyper- 
surfaces in the expanding direction. We further show that the metric g 
decays asymptotically to 7 at a well-defined rate (that correctly predicted 
by linearized theory) and give the sharp rate of decay. In this sense our 
result may be viewed as a nonlinear stability result for the future evolution. 
We could also view it as implying nonlinear instability for the past evolution 
but, since our arguments are insufficient to treat global evolution in the past 
direction, we shall concentrate here on the expanding direction. Since the 
formation of black holes would be expected to violate geodesic completeness 
towards the future, we can also interpret our smallness condition in the data 
as sufficient to exclude the formation of black holes. 

We work in a specific gauge defined by constant mean curvature slicing 
and a spatial harmonic coordinate condition which serves to kill off certain 
second order terms in the spatial Ricci tensor, reducing it to a nonlinear 
elliptic operator on the metric. This in turn effectively reduces the evolu- 
tion equations for the spatial metric to nonlinear wave equations wherein, 
however, the lapse function and shift vector field are determined by an asso- 
ciated set of (linear) elliptic equations. Local existence and well-posedness 
for the Einstein equations in this gauge was established in [T] along with a 
continuation principle which provides the needed criterion for proving global 
existence. 

The main tool we employ for our global existence proof is an energy argu- 
ment based on the Bel-Robinson energy and its higher order generalization, 
which we define. The Bel-Robinson energy for a vacuum spacetime is basi- 
cally an L 2 -norm of spacetime curvature on a given Cauchy hypersurface, 
and its higher order generalization incorporates the L 2 -norm of the spatial 
gradient of this same curvature. One of the key steps in our proof will be 
to show that, in our chosen gauge, this generalized Bel-Robinson energy 
bounds an H 3 x H 2 norm of the perturbed first and second fundamental 
forms of a CMC slice in the spacetime (M,g). 
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Nontrivial spacetime perturbations which preserve flatness are invisible 
to such purely curvature based energies, and this is the reason we have 
been forced to impose an additional rigidity condition upon the hyperbolic 
manifolds that we consider. Already by Mostow rigidity one cannot perturb 
the flat metric 7 = —dp®dp+p 2r y to another flat one by simply deforming the 
hyperbolic metric 7 on M, but there can be more subtle ways of deforming 
7 on M that preserve flatness. These arise whenever (M, 7) admits so- 
called nontrivial traceless Codazzi tensors. Our rigidity requirement is that 
(M, 7) be such as to exclude such tensors — a condition which is known to 
be satisfied for a non-empty set of hyperbolic manifolds. 

The Bel-Robinson energy is of course not a conserved quantity but, to- 
gether with its higher order generalization, can actually be shown to decay in 
the expanding direction for sufficiently small perturbations of a hyperbolic 
cone spacetime. The main source of this decay is the overall expansion of 
the universe which leads to an omnipresent term of good sign, proportional 
to the energy itself, in the time derivative of this energy. A corresponding 
result holds for the generalized energy. The remaining terms in the time 
derivative in general have no clear sign but fortunately can be bounded by 
a power greater than unity of the generalized energy itself. When the initial 
value of the generalized energy is sufficiently small this implies decay to the 
future at an asymptotically well-defined rate and leads to our main result. 

While we shall not pursue this issue here, there seems to be a straight- 
forward way to remove the rigidity constraint and thereby deal with arbi- 
trary hyperbolic M. This involves supplementing the Bel-Robinson energies 
considered here by another non-curvature-based energy called the reduced 
Hamiltonian. As discussed in [Jj this quantity is always monotonically de- 
caying towards the future (even for large data) but bounds at most the rather 
weak H 1 x L 2 norm of the CMC Cauchy data. However this should more 
than suffice to control the finite dimensional space of moduli parameters 
which arises in the case of non-rigid M but is invisible to the Bel-Robinson 
energies. 

Apart from general Lorentzian geometry results such as singularity theo- 
rems and conclusions drawn from the study of explicit solutions, very little 
is known about the global properties of generic 3+1 dimensional Einstein 
spacetimes, with or without matter, and present PDE technology is far from 
being applicable to the study of such global questions, except in the case of 
small data. 

In [3] Christodoulou and Klainerman proved the nonlinear stability of 
3+1 dimensional Minkowski space, i.e., a small data global existence result 
together with precise statements about the asymptotic decay of the metric 
to the Minkowski metric. This proof was based on a bootstrap argument 
using decay estimates for suitably defined Bel-Robinson energies. A central 
element in the proof was the construction of approximate Killing and confor- 
mal Killing fields, which were then used in a way which is analogous to the 
way in which true Killing and conformal Killing fields of Minkowski space 
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are used in the proof of the Klainerman Sobolev inequalities for solutions of 
the wave equation on Minkowski space. 

In still earlier work |Hj Friedrich had proven global existence to the future 
of a Cauchy surface for the development of data sufficiently close to that 
of a hyperboloid in Minkowski space, with asymptotic behavior compatible 
with a regular conformal compactification in the sense of Penrose. This re- 
sult used the fact that the conformal compactification of such spacetimes 
has a regular null boundary (Scri) and exploited Cauchy stability for a con- 
formally regular first order symmetric hyperbolic system of field equations 
deduced from the Einstein equations. Roughly speaking, local existence for 
the conformally regular system can correspond, for sufficiently small data, 
to global existence for the conformally related, physical spacetime. 

Our argument is close in spirit to that of Christodoulou and Klainerman 
but is much simpler than theirs by virture of the universal energy decay 
described above. The source of this decay can easily be seen in linear per- 
turbation theory by exploiting the fact that p-^ is a timelike homothetic 
Killing field in the background. One readily constructs from this an exactly 
conserved quantity for the linearized equations which differs from the (lin- 
earized analogue of the) Bel-Robinson energy we consider by a multiplicative 
factor in the time variable p. This gives immediately the specific decay rate 
predicted by linearized theory and our arguments ultimately show that this 
is the precise decay rate asymptotically realized by solutions to the (small 
data) nonlinear problem. 

Our arguments are also similar in spirit to those of [3] in which Choquet- 
Bruhat and Moncrief treat perturbations of certain £/(l)-symmetric vacuum 
spacetimes onlxEx^ 1 , where £ is a higher genus surface and in which the 
U(l) (Killing) symmetry is imposed along the circular fibers of the product 
bundle RxExS^MxS. Their results also use energy arguments which 
exploit the universal expansion to obtain decay for small data. A significant 
generalization of that work is presented in the article by Choquet-Bruhat in 
the present volume, wherein she removes the restriction to "polarized" solu- 
tions adopted in the earlier work. For the case of linearized perturbations, 
Fischer and Moncrief [HJ have analyzed the stability of higher dimensional 
analogues of the hyperbolic cone spacetimes described above wherein the 
hyperbolic metric 7 is replaced by an arbitrary Einstein metric with neg- 
ative Einstein constant. These of course include the higher dimensional 
hyperbolic metrics but in fact comprise a much larger set. It now seems 
likely that the nonlinear stability problem for these spacetimes can be han- 
dled by a combination of the methods employed herein and in the article by 
Choquet-Bruhat . 

We now give a more precise description of our main results. Let g be a 
Riemannian metric on M and let fcbea symmetric covariant 2-tensor on M. 
We call (M, g, k) a vacuum data set for the Einstein equations if (g, k) satisfy 
the vacuum constraint equations, reviewed in section TTT\ below. Given such 
a vacuum data set there is a unique maximal Cauchy development (M, g) 
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of (M, g, k) which contains the latter as an embedded Cauchy hypersurface. 
Our results concern the structure of (M,g), especially to the future of the 
Cauchy hypersurface, for (g, k) sufficiently close to the data corresponding 
to a rigid hyperbolic cone spacetime (M, 7). We show in this case that, 
in the expanding direction, (M,g) is globally foliated by hypersurfaces of 
constant mean curvature and that (M,g) is causally geodesically complete 
in this (future) direction. In particular, (M,g) is inextendible in the ex- 
panding direction and thus our results support the strong cosmic censorship 
hypothesis. 

Our main result is summarized as follows. 

Theorem 1.1. Let (M, 7) be a compact hyperbolic 3-manifold and assume 
that (M, 7) is rigid (i.e., admits no nontrivial traceless Codazzi tensors). 
Assume that (M,g°,k°) is a CMC vacuum data set with (g°,k°) £ H s x 
H s , s > 3, having to = ti g ok° = constant < 0. Then there is an e > so 
that if 

ll§s°-7llH3 + |||fc - 7 ||„ 2 < e 

then 

(1) The maximal Cauchy development (M,g) of the vacuum data set 
(M,g°,k°) has a global CMC foliation in the expanding direction (to 
the future of M to in CMC time t = tr g k). 

(2) (M,g) is future causally geodesically complete. 

Remark 1.1. (1) Under our conventions, c.f. section® the standard 
hyperboloid {(x,x) = —1} in / + ({0}) C M 3 ' 1 has mean curvature —3 
and Vol(M, g) increases as t /* 0. 

(2) (jfg, yfc) are rescaled Cauchy data that reduce to (7, —7) for the 
background solution. Our energy arguments show that the rescaled 
data approach their background values at a well-defined asymptotic 
rate as t = tr g k /* 0. 

(3) By exploiting the scaling with respect to t at to one can satisfy the 
smallness condition for initial data (g°, k°) corresponding to arbitrar- 
ily large initial spacetime curvature. In this sense one can choose the 
initial hypersurface to be "close to the singularity". 

In outline our paper proceeds as follows. Some preliminaries and a dis- 
cussion of the Einstein equations in our chosen gauge including a review of 
the local existence theorem proven in , are given in sections and 12.11 
Sections I2.2H2.4I discuss the background spacetimes, the constraint set for 
the perturbed spacetimes and the rigidity condition needed to exclude the 
occurrence of a moduli space of flat perturbations. Section 3 introduces 
Weyl fields in the spirit of Christodoulou and Klainerman and presents the 
field equations they satisfy when Einstein's equations are imposed. Section 
0] discusses the Bel-Robinson energy and its higher order generalization and 
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computes the time derivative of these quantities in the chosen gauge. Sec- 
tion describes the scale-free variants of these energies that are used in 
our estimates and section 14.21 gives the calculation which shows how these 
energies actually bound Sobolev norms of the perturbed data in the rigid 
case. Sections El and 15.11 discuss estimates and the differential inequalities 
satisfied by our rescaled Bel-Robinson energies. The global existence proof 
is completed in section El and section 16.11 establishes causal geodesic com- 
pleteness. A number of useful definitions and identities are collected in the 
appendix. 

2. Preliminaries 

Let M be a spacetime, i.e an n+1 dimensional manifold with Lorentz 
metric g of signature — + • • • + and covariant derivative V. We denote by 
(•, •) the scalar product defined by g on TM. Let M C M be a spacelike 
hypersurface with timelike normal T, (T,T) = —1, and let t be a time 
function on a neighborhood of M. Then we can introduce local coordinates 
(t, x l , i = 1, . . . , n) on M so that x l are coordinates on the level sets M% of 
t. We will often drop the subscript t on Mt and associated fields. 

Let dt = d/dt be the coordinate vector field corresponding to t. The 
lapse function N and shift vectorfield X of the foliation {Mt} are defined by 
dt = NT + X. Assume T is future directed so that N > 0. The space-time 
metric g takes the form 

g = -N 2 dt ®dt + g ij (dx i + X*dt) <8> (dx j + X j dt). (2.1) 

Let {ej}j=i v .. jn be a Fermi-propagated orthonormal frame tangent to M t , 
i.e. (V^e^ey) = 0, Vi,j, with dual frame {e*}" =1 . If one drops the as- 
sumption that the frame is Fermi propagated, then in general V^ej = 
e{ + (N^ 1 \7iN)T, where V^/e, denotes the tangential part of V^ej. With 
eo = T, {e^}" o is an ON frame on M, adapted to the foliation {M t }. We 
will use the convention that lower case latin indices run over over 1, . . . , n, 
while greek indices run over 0, . . . , n. Our conventions for curvature as well 
as some useful identities are given in Appendix lA.il The index T in a tensor 
expression denotes contraction with T, for example VxA a = T@V ' pA a . 

The second fundamental form kij of M t is given by k%j = —^{C-T9)ij- In 
terms of the Fermi-propagated frame {e^} we have the following relations 
between N, T and kij . 

= ~ kyT, VjT = —kijej, (2.2a) 

V T e; = (A r_1 VjA / ')T, V T T = (N^ViN)^. (2.2b) 

In computations we frequently make use of equations Q2.2JI to do an n+1 
split, for example VjA,- = VjAj + kijAx- 

2.1. The vacuum Einstein equations. The vacuum Einstein equations 



R a /3 — 0, 



(2.3) 
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consist after a n+1 split of the constraint equations 

R - \k\ 2 + (trk) 2 = 0, (2.4a) 

VjtrA; - V j k i:j = 0, (2.4b) 

and the evolution equations 

£d t 9ij = -ZNkij + Cx9ij, (2.5a) 

£dthj = -ViVjN + N(Rij + tvkhj - 2k im k m j) + Cxhj- (2.5b) 



We will call a solution (g°, kp) to the Einstein vacuum constraint equations 
on M, a vacuum data set. A curve t i— > (g,k,N,X) solving the Ein- 
stein vacuum evolution and constraint equations corresponds to a vacuum 
space-time metric g via (|2.1|) . A vacuum space-time (M,g) with an iso- 
metric imbedding of a vacuum data set (g, k) on M is said to be a vacuum 
extension of (g,k). 

Let g be a fixed C°° Riemann metric on M with Levi-Civita covariant 
derivative V and Christoffel symbol T k y Define the vector field V k by 

V k =g^e k {V i e J -V l e j ). (2.6) 

In terms of a coordinate frame, V k = g^(T k j — f^-). The identity map 
Id : (M, g>) — ► (M, 5) is harmonic exactly when V k = 0, see \Q for discussion. 
A vacuum data set (g, k) is in CMCSH gauge with respect to g if 

tigk = t (Constant Mean Curvature), (2.7a) 

V k = (Spatial Harmonic coordinates), (2. 7b) 

Given a space-time (M,g), a foliation {M t , t S (T_,T + )} in (M,g) is 
called a CMC foliation if VtrA: = for all t € (T_,T + ). In this case, we 
may after a change of time parameter assume t = trk. If the induced data 
(g,k) on M t is in CMCSH gauge for all t G (T_,T+), then {M t } is called a 
CMCSH foliation. The CMCSH gauge conditions imply the following elliptic 
equations for the lapse and shift 

-AN+ \k\ 2 N = 1, (2.8a) 

AX 1 + R l f X-f - C x V l = (-2Nk mn + 2V m X n )e\V m e n - V m e n ) 

+ 2V m Nk i m - V l Nk™, (2.8b) 

where AX 1 = g mn 'S/ m V n X % . The ellipticity constant A[g] of g is defined as 
the least A > 1 so that 

k- x g(Y,Y)<g{Y,Y)<kg{Y,Y), VY G TM. (2.9) 

Let g defined in terms of g,N,X by ()2.1j) . Define A [g] by 

A[g] = A[g] + ||JV|| L « + IIA" 1 )!^ + \\X\\ L oo. (2.10) 

Then 5 is a nondegenerate Lorentz metric, as long as A[g] is bounded. 

We refer to for the background and proof of the following theorem and 
for the analysis concepts used in the present paper. 
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Theorem 2.1 (PQ). Assume that M is of hyperbolic type with hyperbolic 
metric g of unit negative sectional curvature. Let (g°,k°) € H s x IP - x , 
s > n/2 + 1, s integer, be a vacuum data set on M in CMCSH gauge with 
respect to g. Let to = tr/c°. The following holds. 

(1) Existence: There are T^ < to < T + < so that there is a vacuum 
extension (M,g) of (g°,k°), M = (T_,T + ) x M, g £ H S (M), and 
such that the foliation {M t = {t} x M, t E (T_,T+)}, is CMCSH. 

(2) Continuation: Suppose that (T_,T + ) is maximal among all inter- 
vals satisfying point^[ Then either (T_,T+) = (— oo,0) or 

limsup (A[g] + ||Z)^||ioo + ||/c|| L cx)) = oo 

as t y T + or as t~\T-. 

(3) Cauchy stability: Let g be the space-time metric constructed from 
the solution (g, k, N, X) to the Einstein vacuum equations in CM- 
CSH gauge. The map (g°,k°) — > g is continuous H s x H s ~ l — > 
H s ((t_,t + ) x M), for all t-,t+, satisfying T_ < i_ < t+ < T + . 

2.2. Hyperbolic cone space— times. Let (M, 7) be a compact manifold of 
hyperbolic type, of dimension n > 2, with hyperbolic metric 7 of sectional 
curvature —1. The hyperbolic cone space— time (M, 70) with spatial 
section M is the Lorentzian cone over (M, 7), i.e. 

M = (0, 00) x M, 7 = -dp 2 + /3 2 7. 

Let (M,7) be a hyperbolic cone spacetime of dimension n + 1. The family 
of hyperboloids M p given by p =constant has normal 

r = d p . 

Here T is future directed w.r.t. the time function p. Construct an adapted 
ON frame T,ej on M. A calculation gives 

kij = p^J' 

and the mean curvature is given by trk = —n/p. The mean curvature time 
is defined by setting t = trk and the t-foliation has lapse 

N = -(d t ,T) = ^. 

In terms of the mean curvature time we have 

2 

n n 
g(t) = ^7, Ht) = -7. (2.11) 

In the rest of this section we will consider CMCSH foliations, with the 
reference metric g chosen as g = 7. 
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2.3. The constraint set and the slice. Let M be a compact manifold of 
hyperbolic type, of dimension n > 2 with hyperbolic metric 7 of sectional 
curvature —1. 

For s > n/2, let Ai s be the manifold of Riemann metrics of Sobolev class 
H s on M. Then M s is a smooth Hilbert manifold and T> s+l acts on A4 S as 
a Frechet Lie group. 

Lemma 2.2. Let s > n/2 + 1 and fix r G R, r 7^ 0. There is an open 

2 

neighborhood C -M 5 0/ ^2-7, so £/iat /or a// g G W^, i/iere is a unique 
<ft e V S+1 (M), so that (ft : (M,g) -> (M,7) is harmonic. 

Proof. M is compact and 7 has negative sectional curvature. Then there 
is a unique harmonic map (ft £ -ff s+1 (M;M) from (M,g) to (M, 7) 6 . For 
5 close to 7, the implicit function theorem shows (ft is close to the identity 
map Id and hence (ft £ V S+1 (M). □ 

Let be as in Lemma l2~2l Let C 7W S be defined by 

S* = {g e : Id : (M, 5 ) -> (M, 7) is harmonic}. (2.12) 

For g € W*, if is the harmonic map provided by Lemma f2.21 G 5*. 
By uniqueness for harmonic maps with target (M, 7), it follows that 5* is a 
local slice for the action of T> s+l on A4. For s > n/2, let 

= {(g, k) eH s x ff 8-1 , tr 9 A; = r, 

(5, k) solves the constraint equations (|2.4|) }. (2.13) 

be the set of solutions to the vacuum Einstein constraint equations, with 
trk = t. As M is a manifold of hyperbolic type, is a smooth Hilbert 
manifold, cf. 7 j. The action of D s+1 on C£ is the lift of the action on M, 

2 

and therefore the local slice S^. C A4 S lifts to a local slice 3D* , at (^-7, fc) G C®, 

={( 5 ,fc) : ( 5l fc)eq and 5 G <S*}. 

The slice is a smooth Hilbert submanifold of C%. 

A symmetric 2-tensor h on (M,g), which satisfies tih = 0, div/i = 0, is 
called a TT-tensor (w.r.t. g). In the rest of this section, let D denote the 
Frechet derivative in the direction (h,p) G T( 7 _ 7 )C_ n . It is important to 
keep in mind that expressions involving D are evaluated at (7, —7). 

Lemma 2.3. 

T 2 „ S T = {(h,p), h,p TT-tensors w.r.t. 7}. (2-14) 

Proof. We give the proof assuming r = — n, the general case follows by 
scaling. First note 

= Dtrfc = tr 7 /t + ti y p. (2.15) 
Since j = ~7, j235| implies D|A;| 2 = 0. Let # = + (trA:) 2 - |A;| 2 , so 

that the Hamiltonian constraint (|2.4a|) is = H. By the above, 

= ~DH = "DR. 
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Any symmetric 2-tensor can be decomposed as 

h = ipj + /itt + £y7, 

where if} is a function, /itt is a TT-tensor and Y is a vector-field. As 
7 is hyperbolic, R[y] is constant, and due to covariance of R, DR.CyJ = 
YR[^f] = 0. The Frechet derivative of the scalar curvature is the operator 

DR.u = -V fc V fc V + VV 3 '^ - Riju ij , 

which by the above gives 

BR.h = BR.(irf) 

= —{n — 1)A 7 '0 + n(n — 

In view of the fact that A is negative semidefinite, = Di? implies ip = 0. 
Thus, 

h = h T T + C-Yl- (2.16) 
Let V be given by (j2,fij) . Then with 5 = 7, 

J 2 

where V, T, tr are defined w.r.t. 7. By definition, V = on <S r , and therefore 
S r( 7) _ 7 )S r implies using (|2.16j) 

DV = DV.Lyj, 

which by the uniqueness of harmonic maps with target 7 implies that Y = 
and hence h = /itt- In particular tr/i = and therefore by (j2.15|) . trp = 0. 

Let Cj = Vjtr/c — V^kji so that the momentum constraint (|2.4b|) is = 
Cj. By assumption, Vjtr/c = 0, which using h = /itt and the momentum 
constraint (|2.4b|) gives 

= V j p jh 

where V is the covariant derivative w.r.t. 7. By the above, trp = so p is 
a TT-tensor w.r.t. 7. □ 

2.4. Flat space— times. Let (M,7) be a hyperbolic cone space-time of 
dimension n + 1, n > 2, with spatial section (M, 7). Consider a vacuum 
metric 5011M. Then, C a p^s = Rap^s is the Weyl tensor and by the structure 
equations, 

C iTj T = Rij ~ k m k m j + kijtrk (2.17) 

CmTij = d^kmij (2.18) 

where the covariant exterior derivative dFu on symmetric 2-tensors is 

(d v u)ijk = VkUij - VjUik 

Let Eij = CiTjT and F m ij = C m Tij, considered as tensors on M, and define 
the second order elliptic operator A on symmetric 2-tensors by 

Au = V*Vu - nu, (2.19) 
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so that (Au)ij = —V k VkUij — nuij. By ^2 Lemma 4] 

ker A = ker tr n ker d v . (2.20) 

An element of ker is called a Codazzi tensor, i.e. the kernel of A consists 
of the trace-free Codazzi tensors. Clearly, a trace-free Codazzi tensor is also 
a TT-tensor. Let D denote the Frechet derivative in the direction 

(h,p) G T( 7 _ 7 )E T , 

as in section l2~3l 
Lemma 2.4. 

2D£( 7 , -j){h,p) = Ah-(n- 2)h - 2(n - 2)p, (2.21a) 
F>F( 1 ,- 1 )(h,p)=cF(p + h). (2.21b) 

Proof. By Lemma 12.31 are TT-tensors w.r.t. 7. If /i is a TT-tensor, 

w.r.t. 7, then 

Di^-./i = -V*V/t i;? - - nhij = -^Vkhij - nhij, (2.22) 

which gives ()2.21aj) after simplification. The Frechet derivative of is 
given by 

DT) k .h = -g im {Vjh km + V k h jm - V m h jk ). (2.23) 
A computation using (|2.23j) and (|2.18j) yields 

DCtoTY; — Vj/lj m Vj/lj' m + V jPim V ' iPjmi 

which gives 

Consider a curve g\ of vacuum metrics on M, go = 7, such that {Mt} is 
CMCSH foliation with respect to g\, and let g\, k\ be the induced data on 
M_ n . Then 

{h,p) = T7r{g\,k\) 

ax A=o 

satisfy (/i,p) G T( 7) _ 7 \£_ n . As above, let D denote the Frechet derivative in 
the direction (h,p). If we further assume that g\ is a family of flat metrics, 
then DE = and DF = 0. 

Decompose h, p using the L 2 -orthogonal direct sum decomposition ker A® 
ker T ^4, and write h = h° + h 1 , p = p° + p 1 , with hP,p° G kerj4, /i 1 ,^ 1 G 
ker T A Then Di? = is equivalent to the system of equations 

(n - 2)[h° + 2p°] = 0, (2.24) 

Ah 1 - (n - 2)h x - 2(n - 2)p x = 0. (2.25) 

By Lemma f2.3l are TT-tensors on (M, 7). Therefore, by (|2.21b|) . h+p G 
ker tr n ker d v , and hence, by (|2.2U[) . h l +p l = 0. This means that equation 
()2.25j) is equivalent to 

= Ah 1 + (n-2)h 1 . 
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The restriction of A to ker A is positive definite, so it follows that h = 0. 
However, we know that h 1 + p 1 = 0, and hence h 1 = p 1 = 0. Thus we have 
shown that (h,p) = (h°,p°). It remains to make use of 1)2.24(1 . In case n = 2, 
this is trivial, while if n > 3, hP + 2p° = follows. 

By construction, ker T)E n ker DF is precisely the formal tangent space 
T^F(M) at 7, of the space of flat Lorentz metrics F(M) on M. Recalling 
that in dimension 2, TT-tensors are precisely trace-free Codazzi tensors j2j, 
we have proved 

Lemma 2.5. If n = 2, T^F(M) is isomorphic to the direct sum of the space 
of TT-tensors on M with itself, while for n >3, T~W(M) is isomorphic to 
the space of trace-free Codazzi tensors on M . □ 

In case n = 2, M is a Riemann surface of genus > 2, and in this case 
T^F(M) has dimension 12genus(M) — 12, while for n > 3, T^F(M) is trivial 
in case (M , 7) has no non-vanishing trace-free Codazzi tensors, a topological 
condition. This motivates the following definition. 

Definition 2.6. A hyperbolic manifold (M, 7) of dimension 3, is rigid if it 

admits no non-zero Codazzi tensors with vanishing trace. A hyperbolic cone 
space-time (M,j) is called rigid if (M, 7) is rigid. 

A computation |llj shows that (M, 7) is rigid in the sense of Definition 
12.61 if and only if the formal tangent space at 7, of the space of flat con- 
formal structures on M is trivial. Kapovich jlUl Theorem 2] proved the 
existence of compact hyperbolic 3-manifolds which are rigid w.r.t. infinites- 
imal deformations in the space of flat conformal structures. We formulate 
this as 

Proposition 2.7. The class of rigid hyperbolic 3-manifolds (M, 7) (and 
rigid standard space-times (M,^)), in the sense of Definition \2. 6} is non- 
empty. □ 



In this section, and in the rest of the paper, let n = 3. A tracefree 4-tensor 
W with the symmetries of the Riemann tensor is called a Weyl field. We 
define the left and right Hodge duals of W by 



If W is a Weyl field, then *W = W* and W = — * (*W). Define the tensors 
J and J* by 



3. Weyl fields 




(3.1) 



(3.2) 



(3.3a) 
(3.3b) 
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Then 

T* — - T V-v • 



and 

Vr..W..fl„ fl = 

3 



^\pWf5\aP = Z e vw5J* U al3i (3.4a) 



V[ e *W 7 ,s] a( g = -^vimjsJ^ a/3- (3.4b) 

The electric and magnetic parts 2£(W), of the Weyl field W, with 

respect to the foliation M% are defined by 

E(W) a/ 3 = WatfvTT", B(W) ap = *W a ^ v T»T». (3.5) 

The tensors E and B are t-tangent, i.e. E a pT@ = B a pT@ = and tracefree, 
g a PE a p = g a ^B aP = 0. It follows that g^Ey = g" II, ) = 0. 

In case (M,g) is vacuum, i.e. R ap = 0, the Weyl tensor C a ^s of (M,g) 
satisfies C a p^s = R a pj6 the Gauss and Codazzi equations can be written in 
terms of E and B to give 

Vi/cj m — Vjfcj m = e^j B(C)i m , (3.6a) 
Rij ~ kimk™ + kijtik = E(C) tj . (3.6b) 

Note that from the definition (jA. 18j) of , ()3.6aj) is equivalent to cF 7 k m ij = 
—e^Bim. Using the definition ()A.13|) of curl and (|A.19|) we get the alternate 
form of (|3.6aj) . valid if (g, k) satisfies the vacuum constraint equations Q2.4j) . 

-(curlfc)y = B(C)ij. (3.7) 

The following identities relate W,*W,E = E(W),B = B(W), cf. eq. 
(7.2.1), p. 169] 



(3.8) 



WijkT = —(-ij m B m k, *WijkT = £ij m E m k, 

Wijki = —€ij m €kinE mn , *Wiju = —eij m ek£ n B mn . 

The tensors X?tE, Vt-B have the property 

g tj V T E i3 = 0, q^VtBh = 0, 

i.e. the pull-back of VtE, VtB to M is trace-free, but VtE, VtB are not 
t-tangent in general. The following result allows us to express derivatives 
of the Weyl field W in terms of E(W), B(W), J(W), J*(W). See Appendix 
I A. 21 for the definition of div and curl. 

Proposition 3.1. Let E,B be the electric and magnetic parts of a Weyl- 
field W and let J, J* be defined from W by \3.4\j- Then 



divEi = +(k A B)i + J TiT , (3.9a) 
div^ = -(k A E)i + J? iT . (3.9b) 
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V T Eij - curlB - = -N-\VN A S) y - -(£ x k\j + -(trfc)#y - J iT j, 

(3.10a) 

V T flij + curl^y = +iV- 1 (ViV A £)y - -{B x fc) y + -(trfcjBy - J* Tj . 

(3.10b) 

Written in terms of £g t , becomes 
N- l C dt E i:j = +cuTWij - JV -1 (VJV A B)y 

5 2 1 

- -{E x k)ij - ~(E ■ k)gij - -tikEij + N^CxEij - JiTj, 

(3.11a) 

N^Co.Bij = -cuvlEij + iV _1 (ViV A 



b -{Bx k)ij -\{B- k) 9ij - ~(trfc)£tf + N- x CxBn - J* Tj . 

(3.11b) 



Proof. We write V a E a i in two ways. First, 

div£* = V a E ai - N-^NEji. 
Secondly, by (EEJ) and (EOo|) . 

= J TiT + (k A B)i + N-^NEji. 

This gives ()3.9a(l and the argument for (|3.9b|> is similar. To prove (|3.10|) . 
first note the identities 

V k W lTjT = V k Eij - (e u m B mj + e, jX m B mi )k l k , (3.12a) 

V k *W lTjT = VkBij + (e u m E mj + e^Em^k. (3.12b) 

From this we get, after expanding the covariant derivative, and rewriting 
using (j3~H|) 

ei mn VnW mT3 T + e^VnWmTiT = 2(cmlE) ij + 3(5 x fc)y - (trfc)fly, 

(3.13a) 

e^V^WmTjT + e^V^WmTiT = 2(^15)^ - 3(£ x fc)y + (trfc)^-. 

(3.13b) 

The Bianchi equations (|3.4j) imply 

ei mn V T W mni r = 2e i mn V„W r TOT jT - 2J* T , (3.14a) 
e t mn V T *W mnjT = 2e l mn V n *W mTlT + 2J ijT . (3.14b) 
Using (f3~H|) and (|2~2|) we get 

e . m «V T W mnjT = -2V T Bij + 2N-\VN A (3.15a) 
ei mn V T *W mn]T = 2V T E ij + 2N~\VN A (3.15b) 
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Using (|3.15|) . multiplying by |, taking the symmetric parts of (|3.14|) and 
using 1)3. 13|) now gives the identities (|3.10|) . It is straightforward to derive 
(TTTTTt from (l3~Tn|l using (fOTl) . □ 

Given a Weyl field W, the covariant derivative VtW is again a Weyl field. 
Proposition 13.11 gives the following expressions for E(VtW), B(VtW). 

Corollary 3.2. 

E{V T W\j = +curlB - -\(Ex % + -(trfe)£y - J iTi , (3.16a) 



B(y T W)ij = -cmlEij - -{B x % + -(trfc)By - J* Tj , (3.16b) 



-(£x% + -| 
where in the right hand side, E,B,J,J* are defined w.r.t. W. 

Proof. Prom the definition and using 1)2.2(1 we have, taking into account the 
fact that E is i-tangent, 

E(y T W) tj = T^T s T u V u W tljS 

= T v V v E{W)ij - N- x V m NW imjT ~ N- 1 V n NW lTjn 

using (l3~£l) and (TOTT) 

= r"V„£(W% + N-\VN A B{W)) lv 
which gives (|3.16a() using (|3.1(Jaj) . The proof of (|3.16b|) is similar. □ 

4. The Bel-Robinson Energy 

Given a Weyl field W we can associate to it a fully symmetric and traceless 
tensor 

Q(W) a(3jS = W^W// + *W a ^ u *W^ s u . (4.1) 

Q(W) is positive definite in the sense that Q(X, Y, X,Y) > whenever X, Y 
are timelike vectors, with equality only if W vanishes, cf. 4, Prop. 4.2] Let 
E = E(W), B = B{W). The following identities relate Q(W) to E and B. 

Q(W)tttt = \E\ 2 + \B\ 2 , (4.2a) 
Q(W) iTTT = 2(EAB) l , (4.2b) 

Q(W) ijT T = ~(E x E)ij - (fix B)ij + ^{\E\ 2 + \B\ 2 ) 9ij , (4.2c) 

where \E\ 2 = E %3 Eij = \E\ 2 , and similarly for \B\ 2 . From equations (|3.8|) 
and (H^aT) it follows that Q(W) T ttt = if and only if W = 0. The 
divergence of the Bel-Robinson tensor takes the form Prop. 7.1.1] 

V a Q(W) a ^ 6 = W/ 5 u J(W)^ u + W^ u J(W)^ 

+ ^VJWW + *w»»r{w) tl5v , (43) 

and the definition of E(W) and ^(lA 7 ) gives 

V a Q(W) aTT T = 2E l \W)J{W) lTj + 2B ij (W)J*(W) iTj . (4.4) 



16 



L. ANDERSSON AND V. MONCRIEF 



Let W be a Weyl field and let QiW) be the corresponding Bel-Robinson 
tensor. Then working in a foliation Mt, we define the Bel- Robinson energy 
Q(t,W) by 

Q(t,W)= [ Q(W) T TTTdVMf 
By the Gauss law, this has the evolution equation 



d t Q{t, W) = - I NV a Q(W) aTT Tdfi Mt 

JM t 

NQ(W) af3 TT7T a/3 d mt , 



(4.5) 

NQ{W) apTT K af3 d mt , 

M t 

where it denotes the "deformation tensor" of T, 

Kaf3 = V Q T^. (4.6) 

The components of tt in terms of an adapted, Fermi-propagated frame are 
as follows: 

■Kij = -kij, 7r iT = 0, (4.7a) 

7r Ti = JV- 1 V»JV, ittt = 0. (4.7b) 

We will need control of g in H 3 , and for this purpose we consider in ad- 
dition to the Bel-Robinson energy of order zero, Qo(t,W) = Q(t,W), the 
first order Bel-Robinson energy Qi(t,W) = Q(t, VtW"). In the vacuum 
case, J{W) = J*(W) = 0, so we may view Q\ as a function on the set of 
solutions to the Einstein vacuum constraint equations, by using Corollary 
21 to compute E(V T W), B{V T W). 
Expanding d t Q(t,W) using (jO|>. (|OJ) , (|4~7)) . gives 

d t Q(f, W) = -3 / iV[(E x E)-k + (B x B)-k- h\E\ 2 + |£| 2 )trA; 
./Mi 3 

- 2iV~ 1 V i iV( J E A B)i}dfi Mt 



2 I Ml-y.lu, I 



'Mt 

(perform a partial integration and use (|A.15)I ) 

= -3 / N[(E x E)-k+(B x B) ■ k - -{\E\ 2 + |B| 2 )trA; 

JM t 3 
- 2curLE ■ B + 2E- cm\B]dnM t 

-2 [ N(E^JiTj + B^ J* T j)dfJ,M t - (4.8) 

JM t 

It is straightforward to show that this expression agrees with that obtained 
after a direct computation of dtQ(t, W) using (|3.11|) . 

Let r = trA; and specialize to a constant mean curvature foliation {M T } 
in the following. With the discussion in subsection 12.21 as a guide we intro- 
duce the following quantities which vanish when evaluated in the standard 
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foliation on a hyperbolic cone space-time, namely the "trace free" part rr of 



7T, 



Tr a /3 = K a/3 + ^ (g a(3 + T a T p ) (4.9) 



and the "perturbed" part of the lapse, 



N = N - (4.10) 



T 



In the following Lemma we record the form of d T Qi, i = 0, 1 which will 
be used in the global existence proof. 

Lemma 4.1. In a vacuum space time, the Bel-Robinson energies Qj(t, W), 
j = 0,1 satisfy the following conservation laws. 

d T Qo(r, W) = -Qo(t, W) - 3 f NQ{W) aPl5 % a ^T & d^M T 
T Jm t 

+ T NQ(W) T TTTdHM T , 

J M T 



(4.11) 



d T Qi (r, W) = - Qi (t, W) - 2 I Ng 1 {W)dfiM T 

T J M T 

- 3 / NQ(V T W) a ^sfc a ^T^T s d m r (4.12) 

JM t 

5r f 

+ — NQ(V T W) T TTTdfiM T , 
J Mr 



where 



G^W) = E(V T W) ij (J(V T W) iTj + ^£(V T W%) 



+ B{V T wyi{J*(V T W) iTj + ^B(V T W%). 



(4.13) 



In particular, 



J(V T W) iTj + ^E(VtW)h = Tt^V^WaiTj + ^ x E)ij ~^(Bx B) tj 

(4.14a) 

J*(V T W) iTj + ^B{V T W)ij = ^V^Wasrj + 3(E x B\j (4.14b) 
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k rs S7 s (-e ri n B nj ) 

- h rs h K- • -I- h rs h -F, • -I- h rs h m e n e P F, 
tv rv rs i^ i j ~r rv tv sl ±j r j -p rv tv s c r2 raj rip 

- N-^NiVrEij + K s e sl n B nj + k r s e sj n B ni ) (4.15a) 
k rs V s (e ri n E nj ) 

k s k rs Bij + k k s iB r j + k k s e r j ^mA^ n P 

- N-^NiVrBij - k r s e sl n E nj - k r s e sj n E ni ) (4.15b) 

Remark 4.1. In the proof of the main theorem, it is of central importance 
that the terms in \4-14\) are quadratic in 7r,VW, W. This has the conse- 
quence that the term given by j4-13[ ) can be treated as a perturbation term in 
case of small data. In particular the terms fr^V uffl a iTj and fc a ^V ^*W a iTj 
when expanded are seen to be of third order in N~ l ViN,kij,Eij,Bij and of 
second order in kij,V iEj^^V iBj^. We will not make use of the explicit ex- 
pression for dtQi, but for completeness, it is given in equation \4-l^ below. 

Proof. In order to evaluate DwQ(VtW)(T, T, T), we need J(VxW)iTj and 
J*(VrW)iTj- A computation gives 

J{VtW) M5 = V a V T W af3y5 = 7r au V,W af3yS + T v V v J(W)f3 j5 

+ T v R£ a u W^ 5 + T v R^ v W a ^ & (4.16) 
+ T v R^ a u W a ^ 5 + T»R» a u W af3 ^, 

Note that in vacuum, J(W) = and R a fj^8 = W a/ 3 7 s. Substituting R for W 
in (|4,l(i|) and using (|H.8|) to rewrite the terms quadratic in W gives (|4.14a|) . 
A similar calculation for J* (VyVK), taking into account the fact that in this 
case, R = W and *W are distinct, gives 1)4.1 4hj) . It is now straightforward 
to check that (|4.11j) and (|4.12|) hold, given the definition of Q\ in (|4,13|) , □ 

One can decompose equations (|4.14|) into symmetric and antisymmetric 
parts to obtain the analogues of equations (|3.9|) and (|3.1U|) . Setting JiW) = 
J*(W) = for the vacuum case, defining Ey = E{S/xW)ij and By = 
B{VxW)ij and writing E^ and B^ for VkEij and VkBij respectively we 
get, first for the symmetric parts (the analogues of equations ()3.11j0 . 

N^Cg.Eij = k l:j (k lm E lm ) + gij [E lm E lm + k ls k s m E lm + k lm E lm - B lm B lm ] 
+ 2E ij (k lm k lm ) - 3E u E l 3 - 3ki m E mj 
— 3kj E m i — k^kj E[ m + 3BnB j 

+ ~^9jk elkm {Bim\i ~ k ms B s m) + —gik£ lkm {Bj m \i — k ms B s -^ 
+ N- x N\ l E m - h/k s m E im - ^ki'k 8 m E jm 



where 
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+ (trfc)[3^ + h % m E, m + h™E im - 2(trfc)S y - 9ij k lm E lr 



+ N-'N^B^k^ + B r 3 e kir k x k - e rjl B\ - e ril B r 3 \ 



(4.17a) 



N^CqAj = 2B i3 k lm k lm + 9ij [2B lm E lm + k ls k s m B lm + fc' m 5, m ] 



+ N^N^B^ 



— 1 R ' m TP 



3k j Bn ki kj B[ s + kijk B\ m 



~^(kjmk B.U + ki m k B[j) 

— -^9ik^ km {E m j\i — k^Egj^) - -gjk£ lkm (E mi \i — k^ t E si \i) 

+ (trfe)[3Sy + ^fej'By + -jkjBu - 2{tik)B i3 - 9ij k lm B lm ] 

- N~ 1 N\ l [E r i ekj r k k i + E r 3 e kir k\ - e rjl E\ - e nl E r 3 \ 

+ N^CxBij (4.17b) 

and then for the antisymmetric part (the analogues for equations Q3.9J1 ). 



E\ |j = k 3 m B mr ^ r - k 3 m E l Ji + (tvk)e mn k^B sn - k^k^ B SJ e rlj (4.18a) 
B\ lj = -k 3 m E mr e^ r - k^B'J - (tvk)e imn k^E sn + fc* r fc, a V j 

(4.18b) 

The formula for d%Q\ that is analogous to that given above for dtQ is 
given explicitely as follows, 



d f 



d t Qi 



2 j <j Ng jk e 



£km 



k^{B %3 E si \n - E tJ B s 



I'y 



+ N t i g (trk) G _ _„ ( . 

+ fi g N\ £ (E? j E ij \i + & ij B ij \ i 
+ 2Nfi g (k mn k mn ) (tfiEij + B"B, , 

— AN [igkj ^E\El + B\B\ 
+ Nfi g [(E^k, tl )(E mn k mn ) + (B^k^B^k 

- Wn 9 \E i3 E i£ E e 3 - E i:j B ie Bj + 2B i ^E u 
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E)B n e k ^{k k + \g lk (trk)) 



-B)E n e k ^(k M + \g kl {tik)) 



fly 



E tJ Ei m + B l3 B trn 



- 5N» g kfk™(E ij E im + B ij B tm ) | (4.19) 

4.1. The scale— free Bel— Robinson energy. In the rest of section 0J let 
n = 3, and assume (M,g) is a vacuum space-time with a CMC foliation 
{Mr} with r = trfc < 0. Let W a /3 7 s = Ra(3~f6 be the Weyl tensor of (M,g). 
Then W satisfies the homogenous Bianchi identities, i.e. J[W) = J*(W) = 
0. 

Since we will be estimating geometric quantities in terms of Qo, Qi via 
Sobolev inequalities which depend on scale, we need scale-free versions of 
these energies. It follows from the definitions that the following variables are 
scale-free if A has dimensions (length)" 1 . Here indices refer to a coordinate 
frame. 



dab 


\2 

— A 9ab, 




= \- l T a , 


(4.20a) 


kab 




trfc 


= A^tr/c, 


(4.20b) 


N 


= X 2 N, 


Jig 


= aV 9 - 


(4.20c) 



Note that tr/c has dimensions of (length)^ 1 while we treat spatial coordinates 
as dimensionless quantities. The Weyl tensor in the (3,1) form is conformally 
invariant, and hence the Bel-Robinson tensor is also conformally invariant, 
and in particular scale-free. From this can be seen that the Bel-Robinson 
energies Qo, Q\ have dimensions (length)^ 1 and (length) -3 respectively, so 
that the expressions 

Q i = \- 1 - 2i Q i , » = 0,1, 

are scale-free, i.e. Qi is precisely given by Qi evaluated on the scale-free 
variables g, k etc. In the following we will use the scale factor A, defined by 

A = M = _^. (4.21) 
3 3 V ; 

The scale-free energy function which will be used in the proof of global 
existence is the sum of Qo and Q\ , 

S = Qo + Qi- (4.22) 

It is convenient to introduce the logarithmic time a = — ln(— r). Then 
a y oo as r / 0. The logarithmic time a has the property that d a = —rd T 
is scale-free, so that for example d a £ is scale-free. 
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4.2. The Hessian of the Bel Robinson energy. Let M be a compact 3 
dimensional manifold of hyperbolic type and let 7 be the standard hyperbolic 
metric on M. Let (-^7, ^7) be data with mean curvature r for the hyperbolic 
cone space-time (M, 7). For s > 3, let E£ be the local slice for the action 
of T> s+1 on the constraint set C T , at the hyperbolic cone data with mean 
curvature r. 

The energies Qo, Qi,£ may be thought of as functions on the constraint 
set C S T by using equations (|3.6|) and (|3.8|) . The following Lemma is a straight- 
forward consequence of the Sobolev imbedding theorems. 

Lemma 4.2. The scale-free energies Qq(t, W), Qi(t, W),£(t, W), are C°° 
functions on and S^. 

Let HessQo(7> — 7) denote the Hessian of the function Qq, evaluated at 
(7,-7). A computation shows that for (h,p) G T( 7 _ 7 )TS_ n , 

HessQ (7, -7)((M> (M) = \\\Ah\\ 2 L2 + (Ap,p) L 2 + ~\\h + 2p\\ 2 L2 , 

where A is given by (|2.19|) . and || • H^a, (■, *)i a denote the L 2 norm and inner 
product defined with respect to 7. Recall that ker A = {0} if and only if M 
is rigid. It is now straightforward to prove the following Lemma. 

Lemma 4.3. Let M be a compact hyperbolic 3-manifold. The Hessian of the 
scale-free Bel-Robinson energy Qq, defined by equation considered 
as a function on S_3 ; evaluated at the standard data (7, —7), satisfies the 
inequality 

HessQ ( 7 , - 7 )((M, (M) > C(\\h\\h + IWI^O, ( 4 - 2 3) 

if and only if M is rigid. The constant C depends only on the topology of 
M. □ 

Consider a solution h of the linearized Einstein equations on the hyper- 
bolic cone space-time (M,7). The derivative of the Weyl tensor in the 
direction of h, 

W' = DW[y].h, 

is a Weyl field on (M, 7) which satisfies the homogeneous Bianchi equations, 

J(W') = J*(W') = 0, (4.24) 

Let E(W'),B(W') be the electric and magnetic parts of W at M_ 3 C M. 
Then 

DE(V T W)\jh = E(V T W), DB(V T W)\^h = B{V T W'). 

Recall that the second fundamental form of M_3 is k = —7. This implies 
using (EU, (HUH), 

dwE(X7 T W') = &\vB{V T W') = 0, (4.25) 
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which shows that £7(VtW) and B(VtW) are TT-tensors. It follows from 
Corollary 13.21 using k = —7 and (|4.24[) . 

E(V T W') = +cm\B(W') + -txkE(W), 

B(V T W) = -cux\E(W) + -tvkB(W'). 

The scale-free Bel-Robinson energy £ is a smooth function on S^. Using 
the above it is straightforward, using (jA.23() to prove that the Hessian of £ 
is positive definite on H 3 x H 2 in case M is rigid. We state this as 

Theorem 4.4. The hessian Hessf on T, T , evaluated at the standard data 
(-75-7, ^7), satisfies the inequality 

Hess£(£ 7 , ^7)((M, (M) > C(\\h\\m + \\p\\ 2 m ), (4.26) 
if and only if M is rigid. The constant C depend only on the topology of M . 

Results analogous to Theorem 14.41 can easily be proved for even higher 
order Bel-Robinson type energies. This will not be needed in this paper. 

5. Estimates 

In this section we will introduce a "smallness condition" on the vacuum 
data (g, k), under which we are able to control all relevant geometric quanti- 
ties in terms of the energy function £ defined in section 14.11 Recall the def- 
inition of the slice S T in section HOI In particular, vacuum data (g, k) S T, T 
satisfy the CMCSH gauge conditions. 

Definition 5.1. Let (g, k) be a vacuum data set on M with mean curvature 
r. Let A be given by j4-21\ ) and let (g, k) be the rescaled metric and second 
fundamental form as defined in \4.2G\j . Let B{a) be the set of {g, k) £ so 
that 

ll<?-7ll//3 + \\k + l\\ 2 H 2 < a. 
We will say that (g, k) satisfies the smallness condition if [g, k) € B{a). □ 

The smoothness of the scale-free Bel-Robinson energy £ , Lemma 14.21 
together with the fact that the Hessian of £, the scale-free Bel-Robinson 
energy restricted to the slice S r is positive definite with respect to the 
Sobolev norm H s x H 2 , Theorem 14.41 implies, by Taylor's theorem, the 
following estimate. 

Theorem 5.2. Assume that M is rigid. There is an a > so that for 
{g,k) £ B{a), there is a constant D(a) < 00, depending only on a and the 
topology of M , such that 

Diar^ig, k) < \\g - 7^3 + \\k + 7I&2 < D(a)£(g, k). (5.1) 
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In view of the analysis of the elliptic defining equations for N,X in , 
there is a neighborhood of (7, —7) in £^ 3 , such that N, X are small in TV 1 ' 00 
norm, defined by ||/||^i,oo = ||/||l°° + H-D/Hl 00 - It is straightforward to 
check 

Corollary 5.3. Let a > be such that the conclusion of Theorem \5.tH holds. 
There is a constant 5 > so that for (g, k) £ B(a) with £(g, k) < S, it holds 
that 

max(A,\\N\\ w i,oo,\\X\\ w i, aa A\g\\ w i lOB ,\\k\\i,oo) < 1/5. 

Lemma 5.4. Let a > be small enough so that the conclusion of Theorem 
I5.H holds. Let (g, k) £ B(a) and let N, X be the corresponding solutions of 
the defining equations \2. Then there is a constant C such that 

L oo < C£5, (5.2a) 



\\N\\ L °o < C£, (5.2b) 
\\VN\\l°° < C£, (5.2c) 
||#|| L oc < C£s. (5.2d) 



Proof. The inequality ()5.2aj) follows from the definition of B and Sobolev 
imbedding. The Lapse equation ()2.8a|) implies by the maximum principle, 



3 \\h\\ 2 
\N\ < 6 IMl ° 



\k\\ 2 

which gives (|5.2b|) after rescaling. A standard elliptic estimate gives (|5.2cj) . 
Finally, (j4.7j) together with the above estimates yield | |vr| |j^oo < C{£^ +£), 
which after using the smallness assumption and redefining C gives 



\tt\\l°° < C£z . 



□ 



5.1. Differential inequalities for the Rescaled Bel Robinson ener- 
gies. In this section we will estimate the time derivatives d a Qi of the scale- 
free Bel-Robinson energies with respect to the logarithmic time a defined 
in section |4~T1 

Lemma 5.5. Assume (g, k) £ B(a) for a sufficiently small so that the 
conclusion of Theorem \5. 6 A holds. Then 

daQ (a,W) < -{2-2C£ x ' 2 )Q a , (5.3) 

Proof. Replacing all the fields in the RHS of l|4.11|) by their scale-free ver- 
sions, noting in particular that f = —3 and t < 0, we get 



d a Q = -2Q + F 1 
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where F\ is the scalefree version of 

Fi = -9 / NQ abal 7r ab T c T d d mT + 3r f NQtttt 
Jm t jm t 

The maximum principle applied to the Lapse equation (|2.8aj) implies, after 

a rescaling, that N < 3. This gives the estimate 

A<C(||*|| L oc + ||iV|| L oc)Q (r,W). 
To finish the proof note that by Lemma 15, 41 

||£|U» + ll^lk- < c^ 1/2 , 

using the smallness assumption. We write the resulting inequality in the 
form (|5.3|) for convenience. □ 

Remark 5.1. The proof of Lemma \5.5\ gives the inequality 

da Qo < "(2 - C(||#|| L <» + \\ft\\ L o°))Q , 
which is valid without the smallness assumption. □ 

Lemma 5.6. Assume (g, k) € B(a) for a suffiently small so that the con- 
clusion of Theorem \5. 6 A holds. Then 

d a £ < -{2-2C£^)£. (5.4) 



Proof. In view of the smallness condition, the inequality (|5,2b|) and Lemma 
15. 5[ we only need to consider d a Q\ . We proceed as in the proof of Lemma 
15 .5| using the scale-free version of (|4.12|) taking into account f = —3, we 
get 

da-Qi = -2Qi + F 2 
where F2 is the scale-free version of 

F 2 = -6 / Ng 1 (W)d t x Mr ~ 9 / NQ(V T W) abcd 7r ab T c T d d mT 
Jm t Jm t 

+ 5r / NQ{V t W)tttt 

J M T 

We see that in order to estimate F2, we need to estimate Gi(W), which is 
given in Lemma 14.11 Taking into account the detailed structure of Q±, cf. 
Remark 14.11 we get the estimate 

J NGtiW) < c[||^||oo(l + PHooXSq + Qi) 

+ f N\E(V T W)\(\E(W)\ 2 + \B(W)\ 2 )d mT . 

J M T J 

By the Holder inequality, 

\E{V T W)\{\E{W)\ 2 + \B(W)\ 2 ) < CQ\ l2 {\\E(W)\\l A + \\B(W)\\\ A ). 
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By the Sobolev inequality, we may bound the scalefree version of 

\\E(W)\\l 4 + \\B(W)\\l A 
by C£. We now have the estimate 

J iV^! < C ([I^HooCl + ||fc||oo)^ + £" 3/2 ) (5.5) 

< C£ 3/2 , (5.6) 

for (g,k) G B{a). Proceeding similarly with the other terms in F2 yelds an 
inequality which we write in the form (|5,4j) for convenience. □ 

6. Global existence 

Fix To < and let (<?(fo), fc(ro)) be data for Einstein's equations with 
mean curvature To and assume that (g(ro), ^(to)) G &( a ) f° r an a > small 
enough so that the conclusion of Theorem 15.21 holds. 

We have seen above that for small data, the second order scale-free Bel- 
Robinson energy £ satisfies the differential inequality (|5,4[) . We will use this 
to prove 

Theorem 6.1 (Global existence for small data). Assume that M is rigid. 
Let a > be such that the conclusion of Theorem \5.H\ holds. There is an 
e € (0,a) small enough that if (g°,k°) € B(e), then the maximal existence 
interval in mean curvature time r, for the vacuum Einstein equations in 
CMCSH gauge, with data (g°,k°) is of the form (T_,0). In particular, the 
CMCSH vacuum Einstein equations have global existence in the expanding 
direction for initial data in B{e). 
Here e can be chosen as 

e = J D(a)- 1 min(5,C" 2 ), 

where D{a) is defined in Theorem \5. 6 A 5 > is given by Corollarv \5.!A and 
C is the constant in \5.4\) - 

Proof. Under the assumptions of the theorem, by Theorem 15.21 

£{g°,k°) < mm{8,C- 2 ) (6.1) 

holds. Thus the conclusion of Corollary 15 . 31 holds . By the definition of B{a) 
we may apply Theorem 12. II to conclude that there is a nontrivial maximal 
existence interval (T_,T + ) in mean curvature time r, with T + < 0, for 
(g°, k°) in H 3 x H 2 . We will assume T + < and prove that this leads to a 
contradiction, using energy estimates and the continuation principle. 
Let y(o~) be the solution to the initial value problem 

^- = -2y + 2Cy 3 / 2 ,y(a )=y . (6.2) 
da 

Then if yo = £(g°, k°), is such that y(a) < 00 for a G [ao, <r+], we have 

£ (0-) < y(o-), a G [a ,a + ). 
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The solution to (|6.2j) is 

y- 1/2 = C + e^(y 1/2 -C), 

if yo < C~~ 2 , and in this case y(a) < y(cr ) for a G (cro,oo). This means 
that if (|6.1|) holds at cr = (To, it holds for cr G [<ro, <r+). By Theorem 15.21 this 
implies that ||<jf — 7||#3 + + 7||#2 is uniformly bounded for a € [<7o, <r + ). 
By Corollary 15,31 this implies that the inequality 

sup (A[g] + H-DsHloo + PHioc) < <T\ 

ae[(T ,(T + ) 

holds. 

In view of the continuation principle, Point Oof Theorem l2.1| this contra- 
dicts the assumption that (T_ , T + ) is the maximal existence interval in mean 
curvature time r, with T + < 0. It follows that T + = which completes the 
proof. □ 

6.1. Geodesic completeness. 

Theorem 6.2. Let (M,g°,k°) and (M,g) be as in Theorem 16'. it Then 
(M,g) is causally geodesically complete in the expanding direction. 

Proof. By Theorem 16. 1| (M,g) is globally foliated by CMC hypersurfaces 
to the future of (M, g°, k°), i.e. in the expanding direction, with t = trk f 
i* = 0. 

Let c(A) be a future directed causal geodesic, with affine parameter A. 
Let 

u=f y <«,«> = { - 1 , 

be the normalized velocity, where (•,•) = g(-, •). The geodesic equation is 

V u u = 0. (6.3) 
As c is causal, we may use t as parameter. Let 

u° = dt(u) = 

d\ 

In order to prove geodesic completeness in the expanding direction, it is 
sufficient to prove that the solution to the geodesic equation exists for an 
infinite interval of the affine parameter, i.e. limt yt„ X(t) = oo or 

hm / —dt = oo, 

tyt*Jt n dt 



or using the definition of u°, 



/"* 1 

lim / -zdt = oo. (6.4) 

t/t* J t u° 
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Suppose that we are able to prove that Nu° is bounded from above as t /* t*. 
Then (|6.4j) holds precisely when 

lim f Ndt = oo. 

For a function / on M, we have 

df(c(t)) _ d\df(c(t)) _ (dty 1 df( c (t)) _ i 



dt dt dX \dX 



Y L df(c(t)) 1 

) -^x— = ^ Vuf - (6 - 5) 



A calculation in local coordinates using the 3+1 form of g, gives = —N5^ 
where 5^ is the Kronecker delta, i.e. (T,V) = —Ndt(V) for any V. This 
shows that 

-Nu° = (u,T), 

or 

u = Nu°T + Y, (6.6) 

where Y is tangent to M. Let e = 0, 1. Then by assumption, (it, u) = — \e\ 
which using (|6.6|) gives 

\Y\ 2 g = N 2 (u°) 2 -\e\. 

In particular, we get the inequality 

\Y\ g < Nu°. (6.7) 

A computation in a Fermi propagated frame gives using ()6.6|) . 

V U T = Nu N- l ViNei - k tj Y j ei . (6.8) 

By our choice of time orientation we have N > and u° > and tick < 0. 
We now compute using V u u = 0, (|6.6|) and (|6.8|) 

1 



0V2 (u, V n T) 

1 



N(u° 

m (n°VyiV - /,v ; V'V' 



iV(u' 

VyiV , Arf Y i Y'j Ar tr/c 



Nu° t3 Nu°Nu° 3 N 2 (u°Y 



use (|6.7|) and tick < 0, 



< ||ViV|| L co ;fl +[|JVfc|| I/ =c. fl 

use scaling properties of N, k,g, cf. section |4~T1 

< ||VJV||l« ; j + a- 1 I \N~k\\ L °o.~ g , 
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with A = tr/c/3 = t/3. By the proof of the global existence result Theorem 
16.11 we have £{t) < Ct 2 and g is close to 7. Therefore by Sobolev imbedding, 
we can relate the norms w.r.t. g to the norms w.r.t. 7 and we get 

^ln(Nu°) < C (\\VN\\l«> +A- 1 ||iVfc|| L oc) . 

Now an application of the estimate 1)5.2 J) together with the decay of £ gives 

||W||ioo < Ct 2 , 

\\nI\\ lob <ct, 

which in view of \~ 1 t = 3 gives 

ln{Nu°){t) - ln(Nu°){t ) < C, 

and hence ln(Nu°) < C for some constant C as t /* t* = 0. 

We have now proved that Nu° is bounded from above and therefore it is 
sufficient to prove that 

lim f Ndt = 00. (6.9) 

Write N = N + ^ as in PTTO]) . Using (f5~2^) and the scaling rule (jOUj) to 
estimate iV gives N > C/t 2 as t /* t* = 0. This shows that ()6.9|) holds and 
completes the proof of Lemma 16.21 □ 
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Appendix A. Basic definitions and identities 

A.l. Conventions. We begin by recalling some basic facts and definitions. 
We use the following conventions for curvature. 
The Riemann tensor is defined by 

R(X, Y)Z = V X V Y Z - Vy^xZ - V [X ,Y]Z. 

In a coordinate frame {e a } we have 

R d cab Z c = V a V b Z d - V b V a Z d . 

This gives the conventions for index calculations 

[V a ,V b ]t c = R abc d t d , (A.l) 
[V a ,V b ]t c = R c dab t d . (A.2) 
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The Ricci curvature and the scalar curvature are defined (in an ON frame) 

by 

mc(X,Y) = J2(R(ei,X)Y, ei }, (A.3) 

i 

Seal = ^Ric(ej,ej), (A.4) 

i 

or in index notation 

Rij = g kl Rikji, R = g^Rij- 

Note also 

Ric(X,r) = tr(Z i-> R(Z,X)Y). 
The Riemann tensor satisfies the Bianchi identities 

V[eRab]cd = -^{^eRabcd + V a Rbecd + ^ bReacd) = 0. 

The trace free part of the Riemann tensor in an n-dimensional manifold is 

Cabcd = Rabcd TiidacRbd + dbdRac ~ QbcRad ~ dadRbc) 

T (A.5) 

+ 7 7T7 -^{dacQbd — 9ad9bc)R- 

(n — l){n — 2) 

The totally anti-symmetric tensor e in dimension 3 + 1 satisfies the identities 

e aia2Q3Q4 %/w 4 = -det(^;) iJ=1 ,..., 4 , 

fMasat^^ = -det(«5g) ij=2 ,..., 4 , 
e aia2aia4 e aia2 ^ 4 = -2det(^;) JJ=3 ,...,4, (A.6a) 

0:1020304, _ _ef04 

t ta 1 a 2 a 3 f3 i — uu ^ 4 > 

oio 2 "3«4, _ _OA 

c c aia2»3»4 — z, ^ t- 

In an ON frame adapted to a spacelike hypersurface M in a 3+1 dimensional 
manifold we define (c.f. [5J p. 144]) 

f-ijk = f-Tijk- (A. 7) 

e^Siia* = det(<f ) M=1)2 , 3 = e^^Mj 1 , (A.8a) 

^ 2l3 ^ 3 = det(^) fc>J=a> 3) = 2*£<Sg 1 , (A.8b) 

e i ^e hi2h =2S i ^ (A.8c) 

e* 1 ™ 3 *!^ = 6- (A.8d) 
In dimension 3 we have the duality relations 

Cab = £ a b Vm, 
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for Cab = C[ab] ) where 

_ abc 

'Im — 2 fc m Sab- 

A. 2. Operations on symmetric 2-tensors. Define the following opera- 
tions on symmetric 2-tensors on a 3-dimensional Riemann manifold: 

A ■ B = A ab B ab , (A.9) 

(A A B) a = e a bc A b d B dc , (A.10) 

(v A A) ab = e a cd v c A db + e b cd v c A ad , (A.ll) 

(A x B) ab = e a cd e b ef A ce B df + ^(A ■ B)g ab - -(tvA)(trB)g ab , (A.12) 

curL4 a6 = ^-(e a cd V d A cb + e b cd V d A ca ), (A.13) 

divA, = V b A ab . (A.14) 
The operation A is skew symmetric, while x is symmetric, and the identities 
A ■ (v A B) = —2v ■ (A A B) 

A ■ (B x C) = (A x B) ■ C (if tiA = trC = 0) 

hold. The expression Ax B can be expanded as 

(A x B) ab = A a c B cb + A b c B ca 

- -(A ■ B)g ab + ^{tiA){tiB)g ab - {trA)B ab - {txB)A ab . 

A computation shows 

6xv{A A B) = -(curL4) • B + A • (curLB). (A.15) 

Let A be a symmetric covariant 2-tensor on M and suppose A is i-tangent, 
i.e. A a pT@ = 0. Then in a Fermi propagated frame, 

VrAij = TAij (A. 16) 

£d t Aij = NV T Aij -NUkx A)ij + |(fc • A) 9ij 



--{txk){txA)g ij + (trAjhj + (trfc)^ij ) . (A.17) 



2 

"3 ( 

Define the covariant exterior derivative d^u on symmetric 2-tensors by 

(d v u) ijk = V k Uij - VjUik- (A.18) 
The operators curl, div, d v are related by p. 103] 

d V u kij = (cur\u k i + ^(divu, m - V m tru)e™^ i iy (A.19) 

Taking into account the symmetry of curl this implies 

\d v u\ 2 = 2([curlu| 2 + ^|divn - Vtruj 2 ). (A.20) 
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If u has compact support then in dimension n, 

[ \d v u\ 2 = 2[ |Vuj 2 -2 f |divu| 2 + 2 f {u^R k ijl + u ik R k j )u ij . (A.21) 

JM JM JM JM 

This leads to, if tra = 0, 

f (|Vu| 2 + 3Ri jU ki u k J '- lR\u\ 2 ) = [ (|curlu| 2 + ^|dhm| 2 ), (A.22) 
jm 2 J M 2 

in case M is of dimension 3. If we further restrict to (M, 7) with 7 hyper- 
bolic, so that R[y] = —6, we get 

I (|curlu| 2 + ^|dhm| 2 ) = / (|Vu| 2 -3|u| 2 ). (A.23) 
JM 2 J M 
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